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Abstract 

We study the cosmological evolution of an induced gravity model with a scale 
symmetry breaking potential for the scalar field. The radiation to matter tran- 
sition, following inflation and reheating, influences the dynamics of such a field 
through its non minimal coupling. We illustrate how, under certain conditions 
on the potential, such a dynamics can lead to a suitable amount of dark energy 
explaining the present accelerated expansion. 
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1. Introduction 

In induced or spontaneously generated, gravity the gravitational cou- 
pling (constant) and interaction arise as a quantum effect, in particular as a 
one-loop effect in some fundamental interaction, or through spontaneous sym- 
metry breaking, always in association with the coupling of the curvature scalar 
to some hitherto unknown scalar field. Thus gravity itself would not be associ- 
ated with "fundamental physics" but would be an emergent effect in which the 
conventional formulation is the low energy limit. 

Induced gravity models have been applied to cosmology for several years, be- 
ginning with the original model for a time varying gravitational coupling in the 
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presence of matter suggested by Brans and Dicke [2] 0. Subsequently a simple 
scalar field model for induced gravity which avoided the excessive time variation 
of the gravitation coupling was introduced This latter model was globally 
scale invariant (that is did not contain any dimensional parameter) and sponta- 
neous symmetry breaking in such a context not only generated the gravitational 
constant but also a cosmological constant, corresponding to dark energy. The 
introduction of both radiation and matter further showed that the model led 
to Einstein gravity plus a cosmological constant as a stable attractor among 
homogeneous cosmologies and was therefore a viable dark energy model for a 
range of scalar field initial conditions One drawback of the above simple 
model is that it exhibits stability for any (constant) value of the scalar field and 
the desired value must be determined either by the presence of a condensate or 
of quantum effects. 

Subsequently inflation and reheating were studied in an induced gravity model 
having diverse (scale) symmetry breaking potentials @ . In particular the poten- 
tials examined were either associated with the presence of a condensate (Landau- 
Ginzburg) or quantum effects (Coleman- Weinberg). The latter case was studied 
both for an effective potential inspired by the fiat space results and the potential 
obtained in a de Sitter background 0. 

For a potential explicitly exhibiting symmetry breaking, the gravitational con- 
stant is determined by the value of the scalar field at the end of inflation and 
reheating when the field is at the minimum of the potential. At this point, how- 
ever, one no longer appears to have a cosmological constant (or dark energy). 
The scope of this work is to illustrate how a dynamical (time dependent) cos- 
mological constant could arise during the radiation-matter transition. This is of 
interest for the coincidence problem |8[ insofar as the transition from radiation 
to matter will, through the presence of the coupling of the scalar field to the 
Ricci scalar, displace the scalar field from the potential minimum, thus gener- 
ating a cosmological constant. Naturally the cosmological constant generated 
must be such as not to prevent the formation of collapsed structures and the 
return to the minimum be slow enough to respect the current limits on the time 
variation of the gravitational coupling. We shall begin by examining the case 
for a simple Landau-Ginzburg potential and then see how it must be modified 
to fully satisfy the constraints. 

Let us note that besides induced gravity models, there are models wherein the 
term with non-minimal coupling between the scalar field and the scalar curva- 
ture coexists with the standard Einstcin-Hilbert term. In particular, the models 
wherein the inflaton scalar field was non-minimally coupled to gravity have been 
demonstrated to have some advantages with respect to the minimally coupled 
models Q. The hypothesis, identifying a non-minimally coupled inflaton field 
with Higgs boson [10( has allowed to establish relations between the observable 



1 Indeed the technique of the reconstruction of the scalar field potential providing a given 
cosmological evolution, well known for the case of the minimally coupled scalar field, was 
generalized for the case of induced gravity in paper [jj. 
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data coming from cosmology with those from particle physics 11] . While giving, 
in general, more flexibility in tuning the parameters of the model under consider- 
ation, the non-minimally coupled gravity in the presence of the Einstein-Hilbert 
term, loses in comparison with the induced gravity its minimality and an at- 
tractive feature such as a capacity to treat gravity as quantum effect. We shall 
point out later that for our model, the inclusion of the Hilbert-Einstein term 
does not give any advantages and so it makes sense to just consider induced 
gravity. 



2. The model 



Consider the following homogeneous Lagrangian density for a general in- 
duced gravity (IG) model 



1 



(1) 



on a spatially flat RW background ds 2 = dt 2 — a(t) 2 dx ■ dx. 

In the presence of barotropic cosmological fluids, the effective Friedmann and 

Klein-Gordon equations for such a model are 
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(3) 



with the prime denoting a derivative w.r.t. the number of e-folds N = In (a/ao), 
p' M = —3pM (dust) and p' R = —^pr (radiation). We may re- write the equations 
© and ^ in terms o the variables x = ja 2 and y = H 2 as 



V(x) + p M + PR 



and 



y 



y_ 

2y 
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PM 



(4) 



(5) 



I + 67 

where the new variables x and y are positive definite. 

The scalar field a is associated with the observed value of Newton's constant 
and a stabilizing potential V{a) in the form of a symmetry breaking potential 
is generally assumed in order to set its value to ero in such a way that ja 2 , = 
xq = Mp 2 . In the absence of the cosmological fluid, when the fields sits in 
the minimum of the potential, its v.e.v. is different from zero and its zero point 
energy is null. A natural choice for V(a) is that of a Landau-Ginzburg potential: 



V{a) = \ (a 2 



2\2 A 



(6) 
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In a previous paper [B[ we studied an IG dark energy model with just a quartic 
potential (erg = 0). Such a model was shown to be viable as a dark energy 
candidate due to its scale invariance and the presence of the de Sitter attractor 
in the absence of cosmological fluids. On setting er = in ^ the model is 
recovered. These attractors are stable but there is no dynamical mechanism 
selecting among them and fitting observations requires some sort of fine tuning 
of the initial conditions. 



2.1. Absence of Cosmological Fluids 

Let us first consider the case pu — PR = and a generic potential V(a). It 
is clear from (|4I5[) that a de Sitter solution (y — y = const.) exists when 
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The second equation can be re-written as 
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and finally takes the form 
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Eq. (|10p can be solved for x'/x = i 
equation 

12 , 12 7 
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const with 8 satisfying the algebraical 
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The solutions of Eq. (ITT1) are — 
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On substituting such solutions 



into the original equations ([7]) and (JSJ one easily verifies that 62 is associated 
with the potential V oc x oc a 2 . The solution 8\, on the other hand, actually 
corresponds to the trivial case y = and V — (see [3| for details). 
The case 5 = should be studied separately. In this last and Eq. 

© becomes 
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Sx 



dx 
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which is identically satisfied when V oc x oc a , and for more general potentials 
it constrains a;. For the potential ([6]) this last equation leads to: 
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2.2. Matter Domination 

The energy density of the matter fluid enters in the evolution equation of 
the scalar field as an external force which increases the expectation value of the 
scalar field (see Eq. ((SJ)). In the presence of a potential of the form given by 
([6]) one must also take into account the effect of such a potential on the r.h.s. 
of Eq. ([5]) "pulling" the field toward its attractor x n : 



AV - 2x^- = -^-x (x 

ax 7 Z 



xo) => V = 



\xl 



2V 



dV_ 

dx 
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If the scalar field is sitting on the minimum of ([6j at the onset of matter dom- 
ination, its dynamics can be described by two phases. During the first phase 
the field increases because the "pulling" force is negligible at the minimum of 
the potential (V(xq) — 0, dV/dx(xo) — 0) w.r.t. the "pushing" effect due to 
matter, i.e. pu 3> V. One can then approximate Eqs. flU, © by neglecting 
V and its derivatives. During the second phase the energy density of the scalar 
field dominates over matter and the field rolls down the potential and finally 
ends up oscillating with a decaying amplitude around the minimum Xq. If the 
scalar field rolls down slowly enough then it behaves as dark energy and drives 
cosmic acceleration. 

At the onset of matter domination the Friedmann equation 
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well determine the approximate dynamics. The Klein-Gordon equation can be 
cast into the following form 
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where y' jy can be obtained by deriving (1151) and using the continuity equation 
for pM- One obtains 
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and finally Eq. ([TBI) takes the following form: 



The above equation is a first order differential equation for the function x'/x 
and has attractor solutions which can be found by solving Eq. (fT8|) algebraically 
when ^3jQ = 0. One finds 
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The solution — = — 2 was also found in the previous section, in the absence of 
fluids, and still corresponds to the trivial case y — and pu = 0. 
Thus when matter dominates over the energy density of the scalar field, the 
field increases exponentially as 



x = xq exp 



4-< 
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(N - N e ) 
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where N e is the value of N at the matter-radiation equality (we take N = 
today), whereas the matter density decreases as pu = Pm.o ex P ( — 3N) where 
Pm,o is the energy density of matter today. This regime comes to an end when 
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At the beginning of the second phase the above equation can then be recast in 
the following form: 

4 4 

p M ~—V~ —p a (22) 

with e x = (x — xo)/xq, p a the energy density of the scalar field and the 
last equality in (f2"12j) holds approximately only when the scalar field begins to 
dominate and slowly varies in time. In particular if e x < 1 one has p a < pu 
when the scalar field begins to roll toward the minimum of the potential. For 
small values of gamma the condition e x < 1 is generally true and one may 
have a Universe dominated by dark energy today but not in the far past. The 
coincidence problems is then alleviated for such classes of dark energy models 
because reasonable conditions on 7 (small values of 7 are needed to satisfy solar 
system observational limits on the model) put a bound on the maximum value 
of p a . 

On substituting the approximate evolution for x(N) one gets: 



Pm,o , A 
— 5- exp (— 3JV*) ~ r 
Xq 7 2 



exp 



47 
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(23) 



where iV* is the number of e-folds between the epoch when the scalar field be- 
gins to dominate over matter and today. 

When the scalar field dominates, the evolution can be approximated by the so- 
lutions of Eqs. (|4I5|) in the absence of fluids. Such solutions have been obtained 
in Q for the inflationary dynamics and are 



x 1 47 (2 - n) 
x ~ l + 27(n+l) 



(24) 
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and 



y' 47 (2 - n) (n - 1) 
y ~ l + 27(n + l) 



(25) 



where n — . When x'/x <C 1 and y'/y -C 1 expressions (l24l) . (|25|) are very 
accurate and describe a slow-roll dynamics similar to the inflationary case. 

3. Comparison with observations 

IG models as alternative theories of gravity are constrained by many ob- 
servations. The expectation value of the scalar field is associated with the 
"fundamental" Newton's constant Gn, namely the coupling which multiplies 
the Ricci scalar in the classical action. Furthermore the scalar field fluctua- 
tions around its homogeneous/classical trajectory couple to the trace of the 
energy-momentum tensor and they may modify the value of Newton's constant 
measured in Cavendish-like experiments. In particular if these fluctuations are 
massive they mediate short range interactions and the relation x = (8itGn) 1 
still holds at distances larger than their Compton wavelength. On the other 
hand, massless fluctuations of the scalar field, would generate a long range in- 
teraction and correspondently modify the measured Newton's constant by some 
factor 17(7) with lim 7 _>o 3(7) = 1- In the latter case the effective Newton's con- 
stant measured in Cavendish like experiment isCff = (87r;r) 1 5(7). 
The fluctuations of the scalar field around the minimum of the LG potential 
are indeed massive and their mass is m? = 2Xxq. We thus expect modifica- 
tions to General Relativity up to distances of the order r ~ U/Xxq) . Such 
modifications must be compatible with present day constraints on Gjv. New- 
ton's law of gravitation has been extensively tested at millimeter scales, at 
geophysical scales (~ 10 2 m) and at astrophysical scales (~ 10 8 m). The first 
two measurements succeeded in determining the value of the Newton's constant 
today with a precision of about 0.012%. Today we have (87tGjv) _1 = Mp 2 and 
M P ~ 2.436 • 10 18 GeV. 

In IG theories the "fundamental" Newton's constant also depends on time. At 
the cosmological level such a dependence plays a crucial role in the expansion 
of the homogeneous Universe. Given the matter content of the Universe at a 
certain time, larger values of Gn may determine a faster expansion of the Uni- 
verse compared to GR. The Big Bang Nucleosynthesis (BBN) era provides quite 
severe constraints on the expansion rate of the Universe in order to reproduce 
the correct abundances of light elements. Milder constraints are also imposed 
by CMB observations. In this latter case such constraints also depend on other 
cosmological parameters and are more cumbersome to impose. Stringent bounds 
on the time variation of Gn today comes from solar system observations. In 
particular the lunar laser ranging experiments impose the following limit on the 
time variation of Newton's constant: 



H 1 Gn I Gn 



today 



< 0.02. 



(26) 
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Io order to simplify the analysis of the viability of our dark energy model we 
set the minimum of LG potential at xq — Mp 2 . When the scalar field sits 
in such a minimum, deep in the radiation domination era, the value of the 
Newton's constant is that of GR. The bounds imposed on Cm by BBN are 
thus automatically satisfied. We first focus our analysis on the global evolution 
of the Universe from radiation domination until now. We assume a spatially 
flat Universe filled with radiation, matter and the homogeneous scalar field. 
In particular we fix the ratio between radiation and dust energy density today 
according to the best fit of the AC DM model, which will be taken, from here 
on, as our fiducial model. These densities can be independently measured and 
we take 

Pm,o 3520 

in order to reproduce the standard transition between radiation and matter 
domination. This era is crucial in our simulations because at the onset of matter 
domination the scalar field begins to move from the bottom of the potential. 
Given the above assumptions, the free parameters left in our model are 7, A and 
Pm,o- The energy density of the Dark energy component is defined as 

Pde = MIpH 2 - PM - p R (28) 

and today 

p DB ,o ~ 3M 2 H 2 - pm.,0- (29) 

Let us note that, in principle, the definition, given in Eq. (|29[) . does not 
exclude a possibility of having a negative energy density for the dark energy. It 
does not occur in the range of the models parameters considered in the present 
paper. However, the very opportunity of the change the sign of the effective 
density of the dark energy is interesting and and attractive from the point of 
view of some cosmological scenarios, considered, for example, in papers fl2j ]. 

In our simulation we replaced the quantities A and pm.q with /? and A im- 
plicitly defined by 

A = /3A, W f,o = ^-M P , (30) 

where a m is related to the best AC DM estimate of the ratio today between the 
dark energy density and matter density: 



1 Pde 



3 Pm,o 



~ 0.925 (31) 

ACDM 



with 73% of dark energy and 27% of dust. It turns out that the equation 
governing the dynamics of the scalar field is independent of A and one only 
needs to fix 7 and /3 to numerically solve for it. The evolution of the relative 
energy densities Qm , and £Ide, with fii = 3A fl H2 a re independent of A as 
well. We then chose to fix 7 and /3 by comparing the function 
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Pm + Pr + Pde 



»eff = -l + „ PM :* PR (32) 
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Figure 1: Results for the best fit in 7, ft. In the figure on the left we plotted the distance as 
a function of the free parameters 7, 0: larger dots correspond to smaller distances. In the 
figure on the right we plotted w e ff calculated for the best fit 75, /3j, (solid line) and the same 
function in AC DM (dashed line). 



calculated for our model with that evaluated in the ACDM model. In particular 
we minimized the distance between these two functions defined for N e [—12, 0] 
by assuming the following definition of distance between two generic real func- 
tions f(N) and g(N) over the interval N 6 [N t ,N f ]: 



d(f,g) = 



\ 



/£' dN(f(N)-g(N)Y 



Nf - N t 



(33) 



and consequently defining the norm of f(N) as 
comparison gives 7^ = 0.014 and /3& = 0.004 and 



,.,LG „,ACDM\ 
W eff> W eff ) 



,, ACDM I 
"eff 



0.08. 



d(f,0). The numerical 



(34) 



7b and /3 = fib, by min- 
dataset (182 supernovae 



The parameter A can now be estimated, given 7 = 
imizing the % 2 test fitting the supernovae la "gold" 
events). We obtained Ah = 8.3 x 10~ 121 and x m — 165 (# of degrees of freedom 
= 182 — 2 — 1 = 179). In spite of the good agreement with the "w e // test" 
(|38p and supernovae observations the model is quite different from our fiducial 
model. In particular, on assuming a Cosmological Constant energy density of 
Pa = 9.6 x 10~ 121 M P 4 the deviations of H 2 w.r.t. ACDM evolution are close 
to 40% today and the relative abundances of dark energy fluid and dust are 
about 61% and 39% respectively. Let us note that the above value of p\ is that 
which minimizes the x 2 test for supernovae for the ACDM model. On assuming 
the best estimate of p\ {PA° st — 9-04 x 10~ 121 Mp 4 ) the differences between our 
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Figure 2: Results for the best fit in S, /3 for the modified potential with 7 = 5 X 10~ 7 . In the 
figure on the left we plotted the distance as a function of the free parameters <5, /3: larger dots 
correspond to smaller distances. In the figure on the right we plotted the ratio between the 
Newton's constant in General Relativity G^f and that in IG (Gn = (Sivya 2 )^ 1 ) calculated 
with Si,, Pi. 



model and ACDM are smaller but still quite large. 

Serious problems arise if one calculates the value of the "fundamental" New- 
ton's constant today and its variation in time: the maximum displacement of 
the scalar field from the minimum of the potential for 7 = 0(1O~ 2 ) leads to 
a decrease of Gm of about 30%. Moreover its time variation is about one or- 
der of magnitude larger than the present observational bounds. Furthermore 
one can calculate the mass of scalar fluctuations around the minimum of the 
potential and find the Compton wavelength of these fluctuations. It turns out 
that Acompton — 10 23 fcm and thus, for astrophysical scales, Newton's law should 
be modified by taking into account the effect of the fluctuations as well. The 
parameter -yt is then bounded by PPN constraints in modified gravity theories 
and is too large of about 5 orders of magnitude. 

In spite of the several deviations from observations it is worth emphasizing 
that our original, very simple, model with a LG potential is not many orders of 
magnitudes away from experimental data. One may thus expect that with small 
modifications of the shape of the original potential one can actually fit all the 
observations while preserving all the relevant features of the former potential. 
A quite simple modification to the "large field regime" (cr > er ) dynamics could 
be done by multiplying each cto term in the LG potential by some decreasing 
function of M(a) such as M(cto) = 1. On expanding the new modified potential 
Vm(c) around <r one then gets 



o"o + 2a 5(j - ctq ( 1 + 2^-<5ct 



4 



10 



A 2 (* dM\ 2 



(35) 



where is evaluated in a = <7q and is negative. Thus the modified potential 
still has an absolute minimum at <jq and Vm(&o) — as for the LG potential. 
Our aim is to obtain a potential which is closer to Ac 4 in the large field regime 
and we are uninterested in the small field region. 
We chose, for example 

M(cr) = e~^T (36) 

where 8 is a free dimensionless parameter. Let us note that for small and even 
negative a the modification given by (|36l) gives particular shapes to the potential 
which depend on the value of 5. In figure © we plotted the ratio 



where V is either the LG potential or Vm, with M given by (|36|) and different 
choices of 5, and the quartic potential Vq = j<J 4 . The last potential produces a 
de Sitter expansion as the attractor of the dynamics and the ratio (|3T[) , expressed 
as a function of 5a = — — 1, measures the difference between some V and the 
quartic potential itself at a given distance 8a from the minimum cto of V. The 
effect of M is that of decreasing the ratio e(V) for a given displacement 5a. 

The mass of the scalar field fluctuations is proportional to A which is small 
and we thus expect a very light mass. One then needs 7 small enough to satisfy 
the PPN constraints. We thus fixed 7 to its largest value compatible with solar 
system observations (7 = 5 x 1CP 7 ) and let j3, 8 and A vary. Smaller values of 7 
lead to similar results. As in the LG case, /3 and 8 are fixed by minimizing the 
distance (|33|) between w^?P M and w^jj. The parameter A is then determined 
by the supernovae fit. 

We determined /3& = 10~ 12 and Sb = 9.39 x 10~ 7 by first comparing the model 
with ACDM obtaining the corresponding minimum distance 

J/„„M ACDM\ 

{ :Sc»{L ^ SxlO-s. 08) 

IFe// II 

The two functions w^jf and w^j DM are nearly indistinguishable. Note that 
8 is of the same order of magnitude of 7 and this correspondence holds when 
varying 7 for this class of potentials. 

We further find A b = 9.6 x 10~ 121 from SN data corresponding to Xm — 167. 
It is the same value obtained from the supernovae fit in our fiducial model 
confirming the fact that the two cosmological dynamics are extremely close. 
This fact is reinforced by the following results: Qde,o and r2j\/,o are those of 
ACDM model within the experimental errors and the same holds for Newton's 
constant. The time variation of Newton's constant is four orders of magnitude 
smaller then the larger variation compatible with the observational bounds and 
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0.0 0.5 1.0 1.5 

e(V) 

Figure 3: In the figure on the left we plotted the ratio as a function of So for the LG 

potential (solid line) and for the modified potential Vm with M given by (36p with 5 = 1 
(dashed line) and 8 = 10~ 1 (dotted line). The modified potential approaches the quartic 
potential at smaller displacements w.r.t. Vlg- Smaller values of 5 improve the mechanism. 
In the figure on the right we compared the shape of the quartic potential (solid line) to the 
LG potential (dashed line) and the modified potential with 5 = 3. 



the PPN constraints are automatically satisfied by the initial choice of 7. We 
finally note that the Compton wavelength is now Ac om pton — 10 19 km. 
We end by observing that the choice of (155)) is one of the many possible choices 
compatible with the data. The crucial point is that the potential must become 
"scale invariant" as early as possible in order to have a slow-roll so as to minimize 
the time variation of Gn as the field returns to the minimum after the "kick" 
generated by the radiation-matter transition. Such a displacement of the scalar 
field, driven by the matter content of the Universe, alleviates the coincidence 
problem and does not prevent the formation of structures for any 7 compatible 
with PPN observational bounds. It is clear however that in our model the dark 
energy is time dependent and will eventually disappear. 

Let us finally address the possibility of adding an Einstein-Hilbert term to our 
original action (JXJ) with some gravitational coupling Mp such as 

Mp 2 = M|+7crg. (39) 

In this case we have General Relativity plus a non minimally coupled scalar 
field. Due to (f39|) . the coupling between the Ricci scalar and the scalar field 
in this model is smaller w.r.t. ((T|) and when a certain amount of matter is 

added one generally expects a smaller displacement of the scalar field from the 

2 

minimum. In particular for an Higgs- Inflation model such a coupling is ^\ *C 1 
where Co is the v.e.v. of the Higgs and the displacement is tiny for 7 < 1. By 
tuning the self coupling of the LG potential, A, one can set the end of the 
period when the scalar field increases (such a period comes to an end when the 
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"the pushing force" , proportional to R, and the "pulling force" , which depends 
on the self-coupling A, compensate each other). At that time one may have 
a scalar held energy density dominating over matter (and consequently a dark 
energy dominated era) or again a matter dominated era. In the latter case the 
mechanism does not work. This is exactly what happens for the non minimally 
coupled inflation-Higgs models. On the other hand, if one considers the case 

2 

~ 1, then dark energy dominates over matter but no major improvements 
w.r.t. the pure induced gravity case ensues. For this reason, while non minimally 
coupled models generally give more flexibility in tuning the parameters despite 
their loss of minimality, we limited our analysis to the pure induced gravity case. 

Conclusions 

In this paper, in induced gravity, we examined the consequences on cosmo- 
logical evolution of the radiation to matter transition for the case of a scalar 
held at the end of inflation and reheating when such a field is at the minimum of 
the scale symmetry breaking potential. When matter begins to dominate over 
radiation in the early Universe the scalar field is displaced from the minimum 
and begins to contribute to the total energy density which drives the cosmologi- 
cal expansion. In the case of a Landau-Ginzburg potential, for a certain range of 
the parameters 7 and A, the energy density of the scalar field begins to dominate 
over matter at small redshifts thus generating an accelerated expansion. The 
model reproduces a cosmological evolution quite similar to that of the ACDM 
model. For the LG potential, however, the observational bounds on the value 
of Newton's constant observed today are not satisfied. 

Small modifications of the shape of the LG potential in the large field regime, 
in order to make it then more scale invariant, are indeed sufficient to vastly im- 
prove the above result. An example of a modified potential has been proposed 
and studied. Its dynamics is very close to that of ACDM and the model was 
also shown to be compatible with the observational constraints on the current 
value of Newton's constant. 

We finally observe that, in our model, the dynamics of dark energy, triggered 
by the matter content of the Universe, alleviates the coincidence problem in the 
sense that cosmic acceleration is a consequence of, and follows, the radiation to 
matter transition and the energy density driving such an acceleration must have 
a value comparable to the matter energy density. Of course in our approach 
for long enough times the dark energy will disappear as the scalar field slowly 
returns the potential minimum. 
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